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Using an instanton effective action formalism, we compute the critical coupling for the 
nonperturbative formation of a dynamical mass via instantons in non-Abelian gauge theories 
with Nf massless fermions. Only continuous phase transitions are considered. For large values 
of Nf the critical couplings are found to be much smaller than the equivalent critical couplings 
obtained from gauge exchange calculations in the ladder approximation. 
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C3 1 I. INTRODUCTION 



Nonperturbative gauge dynamics present many challenges that are yet to be fully understood. One of these is 
the phenomenon of chiral symmetry breaking. A theory with an SU(N C ) gauge symmetry and Nf fermions in 
the fundamental representation has a global chiral symmetry in the massless limit. Provided that Nf < 5.5N C 
, this theory is asymptotically free. Close to this threshold the coupling is bounded above due to an infrared 
fixed point in the beta function. As Nf becomes smaller the fixed point value increases. Below a critical value 
^\ , for Nf the coupling can grow strong enough to break the chiral symmetry down to the vector symmetry: 

O ■ SU{Nf) R x SU{Nf) L x U(l) v SU(N f ) v x U(l)y (1) 

m ; 

As a result the fermions acquire a dynamical mass, which acts as an order parameter for the spontaneous 
breaking of the chiral symmetry. This much is well understood, but the respective parts played by gauge 
' exchanges and instantons in this process still needs clarification. 

It has been known for some time that gauge exchanges can generate a dynamical mass. To leading order 
(ladder approximation) in Landau gauge the critical coupling associated with gauge exchange dynamics is 

_ 2nN c 
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Assuming a two-loop beta function, one can show that the critical number of flavors, below which dynamical 
chiral symmetry breaking occurs, is p| 



■ N c f =4N c ^ « 4N C . (3) 

f 25N* - 15 W 

It is known that instanton dynamics can also break the chiral symmetry Attempts to investigate 

this possibility have been hampered by the presence of an infrared divergence in the integral over the instanton 
size. In their recent investigation of the critical number of flavors for chiral symmetry breaking via instan- 
tons, Appelquist and Selipsky []5| avoided the problem of infrared divergences by introducing the following two 
modifications of the Carlitz and Creamer gap equation || : 

• They introduced a mass dependence for the fermion determinant that is not only valid for the small mass 
limit H but also for the large mass limit [Q . This function simply makes a sharp transition from the small 
mass behavior to the large mass behavior. The effect is to introduce a fermion mass scale below which 
the fermions are integrated out. 

• In Reference || all couplings in the instanton expression are allowed to run according to the two-loop 
beta function. For those values of Nf where the beta function has a fixed point the coupling remains 
fairly constant until it reaches the fermion mass scale. Fermions are integrated out below the fermion 
mass scale. This causes the beta function to revert back to normal asymptotic free running, which gives 
a suppression of the integrand in the infrared region because of the way it depends on the coupling. 



1 



The result of their investigation was that the critical number of flavors for instantons is about A.77N C . From 
this they concluded that the instanton contribution to chiral symmetry breaking is comparable to that of gauge 
exchanges. 

In this paper we reinvestigate chiral symmetry breaking through instantons, but our investigation differs from 
that of Reference in two essential aspects: 

• The recently proposed instanton effective action formalism ]j| is used to determine the conditions under 
which the symmetric vacuum becomes unstable. This formalism reproduces the Carlitz and Creamer gap 
equation in the small mass limit for SU(2). However, we are not interested in the shape of the mass 
function, but rather in the critical couplings and the critical number of flavors. 

• We concentrate on continuous phase transitions, such as found for chiral symmetry breaking through gauge 
exchanges. The effective action formalism makes it clear that the instanton contribution responsible for 
this type of transition is not the one-instanton amplitude of Figure |l|a but rather the instanton-anti- 
instanton amplitude with two mass insertions, shown in Figure p]b. 



FIG. 1. Leading instanton diagrams: a) one-instanton, b) instanton-anti-instanton with two mass insertions. 

The latter point needs clarification. The effective action can be interpreted as an effective potential. This 
makes it easier to determine when the stability of the symmetric vacuum is affected by the dynamics. For small 
couplings the global minimum of the effective potential is located at the origin, where the order parameter is zero. 
When the coupling is increased the global minimum remains at the origin until the coupling crosses the critical 
value where the phase transition occurs. At this point the global minimum can either jump discontinuously to 
a nonzero value of the order parameter, indicating a first order phase transition, or it can move continuously 
(though non-analytically) to a nonzero value of the order parameter, indicating a continuous phase transition. 
For a first order phase transition a local minimum must develop as the coupling is increased. This becomes 
deeper as the coupling is increased further until it is below the minimum at zero and so becomes the global 
minimum when the coupling crosses the critical value. A continuous phase transition appears when the global 
minimum 'rolls' out of the origin. This implies that the minimum at the origin must flatten out and become 
unstable. To investigate first order phase transitions one must be able to analyze the effective potential at 
arbitrarily large values of the order parameter. Unfortunately the effective potential is not reliable at large 
values of the order parameter because it is not possible to make a reliable truncation of the series of diagrams 
in the effective potential for large order parameters. For a continuous phase transition, on the other hand, the 
order parameter would be arbitrarily small close enough to the transition point. Hence, the order parameter 
can be used as an expansion parameter, allowing one to make reliable truncations of the series of diagrams. For 
this reason we restrict ourselves to the investigation of continuous phase transitions. 

In a stability analysis for continuous phase transitions the one-instanton amplitude is irrelevant. The reason is 
as follows. The effective potential consists of terms with positive powers of the order parameter. Near the origin 
the nondynamical kinetic part of the potential is second order in the order parameter. Hence, to destabilize the 
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potential at the origin the dynamical terms must be of second order or smaller. Terms with higher powers of the 
order parameter will become insignificant compared to the kinetic term when the order parameter approaches 
zero. Looking at the one-instanton diagram, shown in Figure [l|a, one sees that this term always comes with as 
many powers of the order parameter as it takes to close off all the zeromode lines. For Nf flavors there would 
be Nf factors of the order parameter. This means that for Nf > 2 the one-instanton term cannot compete with 
the kinetic term near the origin. The leading contribution to the dynamic part of the potential that is second 
order in the order parameter is the instanton-anti-instanton term with two mass insertions shown in Figure ^b. 
Therefore only this term is considered in our calculation. It turns out that this term does not have an infrared 
problem in our stability analysis. 

The main result of our stability analysis is presented in terms of an expression for arbitrary Nf and N c from 
which one can compute the critical couplings for chiral symmetry breaking via instantons. We determine the 
critical numbers of flavors for N c = 3,4,5 and 6, assuming a two-loop beta function. Numerical values for 
critical couplings are provided in the fixed point regions of these values of N c , up to the respective critical 
numbers of flavors. These results indicate that the critical couplings for the formation of 2-point functions via 
instantons are much smaller than those for gauge exchange dynamics when the number of flavors becomes large. 
The increasing strength of the instanton dynamics can be understood in terms of the increase in the number of 
fermion lines of the instanton vertices. 

The instanton effective action formalism is briefly reviewed in Section |n|. We distinguish between the dy- 
namical and nondynamical terms of the effe ctive action. Section III addresses the nondynamical terms and the 
dynamical term is considered in Section IV. We present and discussed the results of this calculation in Section 
M and conclude with a short summary in Section VI. 



II. INSTANTON EFFECTIVE ACTION FORMALISM 



The 2-point effective action J9|]10| in the instanton formalism of Reference m is given by 

T[S] = Tr {In (S" 1 ) } + Tr { (S' 1 - S' 1 ) S} + W 2PI [S] (4) 

where the fermion modal propagator, S m , consists of different propagators for the different types of modes 
associated with an instanton background configuration. We parameterize the full fermion propagator S by 

SM = 1 = - S fr) (5) 
if — S(p) p 2 + S 2 (p) 

where is the dynamical mass function. 

One can classify the terms of the effective action as either dynamical or nondynamical. Dynamical terms are 
those that vanish when the coupling is set to zero, while the nondynamical terms are those that remain. The 
first two one-loop kinetic terms in (0) are nondynamical terms. They consist of one-loop diagrams that are only 
composed of the full fermion propagator S and the fermion modal propagator, S m . 

The dynamical term W2P1 [S] in (^) is the sum of all 2PI vacuum diagrams constructed with Feynman rules 
for: 

• the 2./Vf-point instanton vertices, V/ and Va, given below; 

• the full fermion propagator S, given in (||); and 

• an integral over the collective coordinates (sizes, positions and color orientations) of all instantons in the 
diagram. 

The dynamics in this analysis are provided by the instantons and are represented by the instanton vertices: 

Nf 

V/ = ^[] (WVS) (fnS^So) , (6) 
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where S 1 is the inverse bare propagator and the functional derivatives Sq(= 5/5t]o) and 6o(= ^/^Vo) operate 
on the zeromode parts of the fermion sources, 770 and rj . The expression for anti-instantons, Va, differs from 

the expression in only in that the zeromode functions, ip° and ip > have the opposite helicities. 
The coefficient of the instanton vertex is the integrand of the 't Hooft instanton amplitude Q : 

A = K1 (a)±r(np) b °+ N f. (7) 
P° 

Here p is a scale parameter for the size of the instanton and 60 = |(lliV c — 2Nf). 
The dependence on the gauge coupling in (Q) is contained in 

2n \ 2N < / 2n 



The factor (pp) b ° appearing in (Q) can be incorporated into the exponent part of 7(a). This then leads to 
one-loop running as a function of p for the coupling in the exponent in 7(a). The instanton expression in (^) 
can be computed to higher orders which would cause the coupling in the monomial part of 7(a) to run as a 
function of p as well, but always at one order less than the coupling in the exponent. To all orders one would 
expect that both couplings run according to the exact beta function. Hence, assuming that the two-loop beta 
function is an accurate enough representation of the exact beta function, one can allow both couplings in 7(a) 
to run according to the two- loop beta function H. This is what we shall assume so that we can make the 
replacement 'y(a(p,))(p,p) b ° — ► 7(a(p)). 

The shape of 7(a) as a function of a is shown in Figure ||. The instanton amplitude peaks at a value of 
a = t:/N c . 




FIG. 2. The shape of 7(a) for N c = 3. 

The prefactor in (Q) is 

2exp(5/6) exp(BN f - CN C }) 

7r2(iV c -l)!(7V c -2)! U 

where we follow || and use the one-loop values: B = 0.2917, C — 1.5114 in MS-bar [{Tit , 

Note that unlike the case in Reference || the integration over the color orientations has not yet been done in 
(Jfy. To alleviate the combinatorical calculations, this step is postponed until after the propagators are connected 
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and unnecessary gauge transformations are cancelled. However, k in (Q) already contains the volume of the 
factor group, SU(N C )/G, where G is the subgroup that leaves the instantons invariant. All that remains of the 
integration over color orientations is a color averaging integral. 

The expression of the fermion zeromode function in singular gauge that appears in (0) is 

*°(«) = , , (10) 

Ky/Jiv X 2 (p 2 + X 2 ) 3 ' 2 

where p is the renormalization scale^ and x 2 denotes an Euclidean dot-product: x^x^. The Euclidean Dirac 
matrices, 7^, are here defined so that {7^,7,,} = 25^ ■ 1, tr{7 M 7„} = 4<5 M „ and 7^ = 7^. The U a denotes a 
Dirac spinor with isospin index a. (The isospin index is implicit in ip°(x).) This spinor is normalized such that 
E a U a U a = 1 and it obeys the identity £ a U a U a = |(1 + 75). 

We shall also later need the Fourier-transform of the zeromode function 



*°(p)= / ^(s)e«p(fe.p) d^ = i ^PPW f IPyUL), (ll) 



;2ir V2p p^ U a t ( p^/p 2 
where the function, 

f{x) = 2x [Io^K^x) ~ h{x)K {x)] - 2I 1 (x)K 1 (x), (12) 
is defined in terms of modified Bessel functions iQ . 

III. THE NONDYNAMICAL TERMS 

The trace over the inverse modal propagator in (Q) contains an integral over all the collective coordinates 
involved in the diagram. The general form of the result of this integration is the inverse bare propagator, 
Sq 1 = if 1 , multiplied by a momentum-dependent scalar form factor, which we shall assume to be equal to 1. 

In terms of (j^) the nondynamical terms in (^), which we denote by Tkim becomes 

r - = -/(?fw- 21 " [!,2 + E2(p,1 )(0- (13) 

It is customary to introduce an infrared cutoff at the momentum scale where po = Yi(po). Close to the phase 
transition point £(p) is very small so that po w 0. An expansion of Tkin in term of T,(p)/p contains only even 
powers. We drop the terms that are independent of S(p) because they do not affect the location of the minimum 
and we drop the terms that are higher than second order in S(p) because they are suppressed. The result is 



4i=-^y o ^i 2 /-/,. (14, 



Now we make a redefinition of the integration variables and do a Fourier transformation. For this purpose 
we define 

/oo 
<j(u>) exp(iivt) du> (15) 
-00 

where p is the renormalization scale. Now we have 

jij jij 4 poo poo poo 

= T^2 — / / < J ( UJ i) (J ( UJ 2) / exp(iujit) exp(iLu 2 t) dt duji duj2 

^71 J —OO J —OQ J — OO 

N c Nfp 4 [°° 

' <7{uj)<j(—uj) auj. (16) 



2n 



x The p is included here to get a mass dimension of | for the fermion zeromodes. 
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IV. THE DYNAMICAL TERM 



The instanton-anti-instanton term that is quadratic in the dynamical mass looks as follows: each mass inser- 
tion closes off two lines of an instanton vertex. The remaining lines are all connected between the two vertices. 
(See Figure [j]b.) The resulting diagram is expressed as 

= Tk 2 f f 7 (a(pi)) 7(«(pa)) ( P 2 P i P2 ) N ? D( Pl )D(p 2 )G( Pl , P2)XX d Pl dp 2 (17) 

J J Pi P 2 

where P i and P2 denote the sizes of the two instantons; the combinatorial factor T simply counts the number 
of terms associated with this diagram; D( Pl ) and D( P2 ) are the parts of the expression due to the two mass 
loops and G(pi, f> 2 ) is the part of the expression due to the remaining fermion loops. 



A. The combinatorics and color averaging 

To count the number of terms we start by considering all the possible ways that the propagators can close 
off the lines of the two instanton vertices to give a diagram with two mass loops. Then we perform the color 
averaging and add up all the different terms. 

Recall that each line of an instanton vertex has a specific flavor associated with it. There is one incoming and 
one outgoing line for each flavor. Because the propagator is flavor diagonal it must connect lines of the same 
flavor. Hence, there are Nf distinct ways to connect one of the two mass insertions, each leaving the other mass 
insertion and the remaining fermion lines with only one way to be connected. 

Next we perform the averaging over the color orientations of the two respective instantons. The color orien- 
tations are represented by gauge transformations that operate on the zeromodes: 

V/^gV and if -^Ijfg' 1 . (18) 

In the mass loops they cancel because there we have a g and a g^ 1 from the same instanton. The gauge 
transformations of the remaining lines can be combined into relative gauge transformations of one instanton 
with respect to the other: g^ x g 2 — g r or g 2 l g\ — g^ 1 - Thus the two color averaging integrals become one trivial 
(= 1) and one nontrivial averaging integral. The nontrivial integral contains (Nf — 1) factors of g r and g^ ■ 
The result of the averaging is a series of tensor structures each multiplied by a coefficient that depends on Nf 
and N c . General expressions for these coefficients are provided in Reference Jl3| . The tensor structures indicate 
how color indices of the same instanton are interconnected. This in turn determines how the fermion lines 
between the two instantons are closed off into fermion loops and thereby determines the form of the momentum 
integrals. 

The evaluation of all these momentum integrals is too challenging to attempt. However, we expect that 
the sizes of the momentum integrals do not differ by more that an 0(1) factor and that t he va riations would 
average out. Therefore we set them all equal and estimate their typical size in Section [V C[ Each tensor 



structure consists of (Nf — 1)! terms because of all the possible permutations of the (Nf — 1) pairs of fermion 
lines between the two instantons. The different tensor structures denote different relative permutations of one 
line with respect to the other in the (Nf — 1) pairs of fermion lines. All tensor structures that are associated 
with the same conjugacy class in the permutation group have the same coefficient. 

Because we assume that the momentum integrals are the same, we only have to add all the coefficients of all 
the tensor structures multiplied by the number of terms associated with each tensor structure. It can be shown 
that for permutations of (Nf — 1) elements, 



Cl 



(N c + N f - 2) 



where Ni is the number of group elements in the i-th conjugacy class of the permutation group and CKj is the 
coefficient of the tensor structures in that class. To get the total count we have to multiply (|l9| ) by Nf for the 
number of ways to connect the mass loops and by (Nf — 1)1 for the number of terms in each tensor structure. 
The resulting combinatorial factor is 



(i 



jV/!(^-l)! 

(JV c + JV>-2)!' ( ' U) 



B. The mass insertions 
Using (||) and we obtain the following expressions for the mass loops 

D(p) = / ^ (p)So 1 (p)S(p)S Q 1 (p)^(p)- d " P 



(27T) 

P 2 /" S 0?) ,2 /PP" 



/ 2 £ *»* 



/i y p 2 + s 2 (p) J v 2 

S(p)/ 2 (f )pdp + 0(S 3 b)/p 3 ) (21) 



P 2 



p 



where we again made use of the approximations introduced in Section [II. Next we again make a redefinition 
of the integration variables and do a Fourier transformation as in ( jl5j ) . In addition we define 



As a result we have 



^=exp(77) and / (exp(a)) = / Q (a) (22) 



/OO pOO 
/ ct(cj) exp(iwi)/Q (t + rj) exp(i) dt duo 
-oc J — OC 

/•OO 

= 4exp(?7) / ct(uj) exp(— iu)r])F(uj) duo (23) 



where 



F(cj) = / exp(iwa)/, 2 (a) exp(a) da. (24) 



C. The remaining fermion loops 



The rest of the fermion lines between the two instanton vertices combine into 

f 4 _ , / 

G( Pl ,P2)= II [y ^*°(p^ - 1 ( Pi )S'(p i )^- 1 fe)* ^)J (2^ 4 [ i> 



(25) 



The fermion lines pair up into fermion loops due to the color contractions that appear after the averaging over 
color orientations. Considering one such loop, we find 

Z(p 1 ,p 2) ^g)=-^(p)5 -Hp)5(p)5 - 1 (p)^(p)^(g)V(9)%)^ 1 (5)*°(^ 

2(2tt) 4 2 2 P-Q f fPlP\ f fP2P\ f (Pl1\ tfP2l\ , n/V 2 / 2\ / or n 

where in the first line we explicitly show the color indices on the zeromodes. By defining 

f d A p 
J(k,p!,p 2 )= J I(pi,p 2 ,p,p-k)-p—s±, (27) 
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one can write (Ea) as 



G{pi,P2) 



Nf-l 

n 



J(ki,pi,p 2 ) 



(2^) 4 <5 4 



(28) 



We shall not attempt to evaluate (]28|) exactly, but instead estimate its size in terms of the number of fermion 
loops. It turns out that the p-integrals are dominated in the region where p\ w p 2 . Therefore we set p\ = p 2 = p 
in the arguments of the /-functions. Then (^) becomes 



P 



(29) 



For the region where p > k we set \J (p — k) 2 = p inside the /-function and for p < k we set v/(p~— k) 2 = k. 
Upon evaluating the angular integral we arrive at 



J{k,pi,P2) 



{^?p\pl 



o 2fc 2 



p dp - 



*1 
2p 2 



/•(f),* 



(30) 



If we assume that p x > p 2 then p — pi, because the /-function with the largest p dominates. Defining 

PPi , fcpi 

x = and y = 

2 y 2 



(31) 



one can write 



J(k,p!,p 2 ) 



P 2 



f 2 {x) x 3 dx + 16 / f 4 (x) x dx 



-f\x) dx 



^P\ 



R(y). (32) 



y 2 Jo ' 

The size of the area under R(y) is ~ This function determines the momentum dependence under the 
remaining momentum integrals in fl28|). To assess the size of the remaining integrals we note that ( p8| ) has 
(Nf — 2) integrals over ki. This will give (Nf — 3) angular integrals. (One can arrange it that each ki at most 
contracts with its two neighbors) . Together with the definitions in ( pl| ) we estimate the size of the remaining 
momentum integrals to be 



N f -1 

n 



d 4 h 



-R 



hpi 



(2^) 4 <5 4 h 



Nf-l 



-i (2tt)- 



Nt—2 



(33) 



Our estimate for the sizes of G(pi,p 2 ) then gives 

G(pi,p 2 ) ~ 



P 4 tt 2 



2pj 

2 4 

M Pi 



(34) 



V. RESULTS AND DISCUSSION 



The integrand of ( |l7| ) is invariant under an interchange of p\ and Therefore, one can split the ^-integration 
into two regions, divided by the value of p\. First, we consider only the part where p 2 < pi. When we substitute 
(B4T) into this part we find 



< 2 < P1 = ^^«V / / ^P 2 ^ D( Pl )D( P2 )p^ 3Nf pf dm ,1,,,. (>',) 



One can see that the ^-integral is dominated by the region near the pi-boundary. For this reason one can 
set j(a(p2)) = 7(a(pi)). The integrand of the pi-integral has its dominant region around some value p m , close 
to \i. Since we consider cases with large numbers of flavors, the coupling runs slowly. Therefore one can set 

7(a(/9i)) = 7(<*Cfiw)) = 7o- 

Next we use ( |22|) to redefine the variables in (jS^) and we substitute in (£3|). The result is 



W, 



(2) _ 



P2<pl 



-j poo roc 

- M 4 ^ K 2 7T 2 2^ 7o 2 / / <7(wi)tT(w 2 )J , (wi)f , (w 2 ) 

*/ — GO J — OO 

/OO /.7/1 
/ exp([5 — 3Nf -j- iuJi\r]i) exp([37V/ — 5 + zu^]^) dr}2 drji duj\ duj2 
-oo J —OO 



(36) 



-oo — OO 



Now we add the part where p 2 > p\. This is the same as 
evaluate one of the w-integrals. The result is 



3N f - 5 + iuj 2 
but with u>i and 0J2 interchanged. Then we 



W, 



(2) 



IP I 



w. 



(2) 



P2<Pl 



w, 



(2) 



pl>P\ 



(3N f - 5)F(c)F(-c) 
(37V/ - 5) 2 +ui 2 



dcu. (37) 
0.3. We notice that 



The maximum value of the part of the integrand in the brackets is at to = with F(0) 2 
this term does not have an infrared divergence. 

Comparing the result in ( ^7|) with the nondynamical term in ([l6]) one can see that the potential becomes 
unstable at the origin when 



^ 2 7T 3 2 Ar / 7o 2 



This can be expressed as 



where 



S(N c ,N f ) 



F(0) 2 
3N f -5 

£70 > 1, 

2F(0) exp(5/6) exp(BN f - CN C 
(N c -2)\ 



> 



N c N f 
2ir 



2 N f+ 1 (N f - 1)! 



(3AT / -5)(iV / + Ar c -2)!iV c ! 



(38) 



(39) 



(40) 



only depends on N c and Nt. 

The critical coupling for instanton dynamics is the one for which 

1 

7 (a(p m ) = 



S(N c ,Nf) 



(41) 



It necessarily has to be smaller than n/N c . The requirement in (f4l|), together with (^0|), is the main result of 
our analysis. From it one can compute a], for any value of N c and Nf > 4/3. Numerical estimates of ct\ for 
various values of iV c and Nf are provided in Figure ^. 

Before we discuss these critical couplings we need to clarify the ranges of numbers of flavors for which these 
critical couplings are computed. These ranges are defined under the assumption that the running of the coupling 
is governed by the two-loop beta function - an assumption which becomes more accurate as the coupling becomes 
smaller. We only compute the critical couplings in the region of Nf where the two-loop beta function has an 
infrared fixed point. The left-hand end of each curve in Figure ^ is therefore the point where a fixed point value 
appears, which moves down from infinity as a function of increasing Nf. 

The right-hand ends of the curves in Figure [| indicate the critical numbers of flavors Nf above which chiral 
symmetry breaking does not occur. When there is a nontrivial infrared fixed point in the beta function the 
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maximum coupling that can be reached is the fixed point coupling. Assuming a two-loop beta function the 
fixed point coupling is given by 



a - b -MnN c -2N f )N c 

c 34/V| - l3N?N f + 3N f y ' 

where b and c are respectively the first and second coefficients in the beta function. The critical numbers of 
flavors Nf is reached when the fixed point coupling given by ( ff2|) falls below the critical coupling, obtained 
from solving ( f4l| ) . Therefore for values of Nf above N* the coupling is unable to reach the critical value so that 
chiral symmetry breaking does not occur. The critical number of flavors, Nj are 13.48, 17.95, 22.41 and 26.87 
for N c = 3,4,5, and 6, respectively. This gives Nj i=s A.5N C , compared to 4.77 N Cl reported in Reference ||. 




FIG. 3. The critical couplings for N c — 3, 4, 5 and 6, each computed over its fixed point region up to its critical number 
of flavors. 

We computed critical couplings for iV c = 3,4,5 and 6 and for Nf in the fixed point region up to Nj. The 
results are shown in Figure^. From these results we make the following observations: 

• The critical coupling decreases as either N c or Nf increases. 

• The values of the critical couplings near Nj are small compared to the equivalent critical coupling for 
gauge exchanges given in (||). Compare, for example, the N c — 3 values. We obtain a ! c = 0.37, while the 
equivalent gauge exchange value is a c { e — ir/4 = 0.785. 

• Toward smaller values of Nf the instanton critical couplings become more comparable to their gauge 
exchange counterparts. 

On this basis we conclude that instanton dynamics are much stronger than gauge dynamics for large numbers 
of flavors. 



10 



VI. SUMMARY 



We have used the instanton effective action formalism of Reference M to compute the critical couplings for 
the formation of 2-point functions. For this analysis we restricted ourselves to continuous phase transitions 
for chiral symmetry breaking, which implies that we only consider the instanton anti-instanton amplitude with 
two mass insertions. We found that this amplitude does not suffer from an infrared divergence. Couplings are 
allowed to run according to the two-loop beta functions. The resulting critical couplings indicate that instanton 
dynamics are much stronger than gauge dynamics for large numbers of flavors. 
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